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Abstract: We first devise a scheme to perform a universal entangling gate via controlled 
collisions between pairs of atomic qubits trapped with optical tweezers. Second, we present 
a modification to this scheme to allow the preparation of atomic Bell pairs via selective 
excitation, suitable for quantum information processing applications that do not require 
universality. Both these schemes are enabled by the inherent symmetries of identical 
composite particles, as originally proposed by Hayes et al. Our scheme provides a technique 
for producing weighted graph states, entangled resources for quantum communication, and a 
promising approach to performing a "loophole free" Bell test in a single laboratory. 

Resume : 

[Traduit par la redaction] 

1. Introduction 

Entanglement plays an indispensable role in many quantum information processing tasks, such as 
long-distance quantum communication [1], teleportation-based quantum computation [2,3], and one- 
way quantum computation [4]. While great progress has been made entangling arrays of neutral atoms 
in optical lattices en masse [5], the current approach to generating such massive entangled states (via 
cold collisions) necessitates state-dependent traps [6]. This state-dependency results in increased noise 
sensitivity and decoherence of atomic qubits [5]. Other proposed approaches for entangling neutral 
atoms feature encodings in vibrational rather than internal electronic states of atoms [7, 8], but are 
subject to similar dephasing of qubits. Approaches based on atomic interactions other than ground 
state collisions have been suggested [9, 10], but none have been successfully implemented and atomic 
collisions still hold the most promise. Thus, there is a need for coUisional quantum gates that allow 
more flexible encodings in robust electronic states — such as the clock states of Rb, Cs, or Group II 
atoms — that are held in state-insensitive traps to minimize decoherence. 

In this work, we examine schemes to entangle pairs of bosonic atoms, analogous to the recently 
proposed fermionic spin-exchange gate [11]. Gates based on this exchange interaction offer a natural 
resistance to errors and more flexibility due to inherent symmetrization conditions. Furthermore, this 
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exchange interaction allows the design of entangling operations for atoms with state-independent (e.g., 
Rb [5]) or partially unknown interaction strengths (e.g., Yb [12] or Sr [13]). The underlying exchange 
interaction for these gates has recently been experimentally demonstrated using bosonic Rb atoms in 
a double- well optical lattice [14]. However, a verifiable entangling gate between an individual pair of 
trapped neutral atoms has not yet been demonstrated. Here, we provide a detailed analysis of these 
operations as they may be carried out using a pair of individually controlled atomic qubits trapped via 
optical tweezers. 

Our approach builds on disparate proposals and experiments for preparing individual atoms from 
a Bose-Einstein condensate [15-17], encoding qubits into long-lived electronic states, coherently ma- 
nipulating and transporting atoms using optical tweezers [18, 19], and performing two-qubit operations 
on pairs of atoms via collisional interactions [11, 14]. The combination of these elements allow for the 
design of a tunable two-qubit gate, which can create an arbitrary degree of entanglement between a 
pair of atoms. We also examine a scheme that exploits symmetrization rules to produce Bell pairs via 
selective excitation. 

These entangling schemes may be realized using qubits stored in the electronic states of a pair of 
atoms trapped with moveable optical tweezers. Trapping at a "magic wavelength" makes the light shift 
potential state-independent. Encoding in atomic clock states — which are insensitive to fluctuations in 
the trapping field — avoids dephasing and ensures qubit coherence during the transport process. Unlike 
the case of a state-dependent optical lattice in which it is trivial to separate the atoms after interaction, 
we have state-independent potentials in which the system's dynamics determine the likelihood of the 
atoms being separated into opposite wells. Under adiabatic conditions, atom separation is guaranteed. 
We consider only the 1-D case for simplicity. Multidimensional effects such as trap-induced resonances 
cannot be captured by the 1-D delta-potential employed here [20] but could potentially be used to 
enhance the atomic interaction further. 

2. Hamiltonian for identical particles in separated tweezers 

The Hamiltonian for two atoms with internal structure in a pair of optical dipole traps (a.k.a., 
"tweezers") is given by, 

(22 \ 

^ = E 1 + "^^^^^ ^(^b' + 2a,,-;k^^^(xa - Xb) U (1) 

2,J=0,1 I ^ ^ J 

where and x\y are the positions of atoms a and b respectively, and are similarly the momenta, 
aij is the state-dependent scattering length that depends on internal atomic states |i)a and |j)b (using 
= |i)a0|j)b)5 ^± is the harmonic oscillation frequency due to transverse confinement [21], and 
d is the time-dependent centre-to-centre distance between wells. For Yb, Sr, and alkali atoms, one 
can usually choose a particular trap-laser wavelength (the "magic wavelength") so that the light shift 
potential becomes state-independent and each atom sees a double-well potential: 

V{x,d) = -y„e-(--^)'/2-' - Ke-(-+^)'/2-'. (2) 

Here, > is the depth of each Gaussian well and is the variance. 

The first three vibrational eigenstates of a single particle in this double-well potential are shown 
in Fig. 1 for varying d. In general, the single-particle eigenstates are {\il)^{d))^ \ilj^{d))^ \ilj^{d))^ . . .} 
and (i-dependence is assumed implicit (e.g., |?/^^) = \ilj^{d))) for notational simplicity. Note that as 
d increases, \ip"^) and become spatially delocalized and energetically degenerate. Thus, when 
\d\ (7 we can write | V^^) = (1^^^) — |^^))/ to represent a single particle localized in the ground 
state of the left well, and similarly = (|V^^) + |V^^))/a/2 for the right well. 

When a second particle is added to the double-well potential, the interaction term in the Hamilto- 
nian may be treated as a perturbation. Accordingly, the new two-particle eigenstates may be written as 
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Fig. 1. The first three eigenstates of a single particle in a double-well potential for different well separations ( 



a sum of perturbed tensor products of one-particle states. We will use a tilde to denote the perturbation 
to the terms composing the new symmetrized eigenstates. For a repulsive interaction between atoms 
(aij > 0), the first six two-particle eigenstates are (see Fig. 2a), 



d = 



d > (T 



(3a) 
(3b) 

(3c) 

(3d) 
(3e) 

(3f) 



States that are antisymmetric under exchange are not affected by the interaction at any separation 
and the tildes have been intentionally omitted from these states. States with atoms in opposite traps 
(e.g., 3a, 3f) are obviously not affected by the interaction in the limit d ^ oo. Note that there is the 
usual on-site interaction penalty for putting two atoms in same trap, resulting in an energy splitting at 
d » (T between states having atoms in opposite traps (3e, 3f) and those having atoms in the same trap 
(3c, 3d), as shown in Fig. 2a. 

In case of attractive interaction (ay < 0), the eigenstates are (see Fig. 2b), 



d = 



d :^ a 

i(I^^^V^) + l^^c^) + I^Js^^ + 1^3^)) (4b) 

^{\r^R) + IV^)) (4c) 

^(|^/.^^/.«) - IV^^/.^)) (4d) 

_ 1^;^)) (4e) 

^(1^^^^) + I W*))- (4f) 
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Fig. 2. Adiabatic energy levels as a function of well separation d for (a) aij — 0.1a and (b) aij — —0.1 a. Well 
separation is in units of a. Energies are in units of hcuo, where cOo is the harmonic oscillation frequency of one 
atom in the ground state of a single well. Symmetric vibrational eigenstates are shown in blue, antisymmetric in 
red. Eigenstates for (a) and (b) are given by Eqs. (3) and (4), respectively. Notice that crossings between oppositely 
symmetrized states are unavoided because the Hamiltonian is symmetric. 



The order of states in (4) is different from (3), since states with atoms in the same trap now posses 
a lower energy due to the attractive interaction. 

Until this point we have neglected the internal structure of the particles (i.e., the qubits). The 
eigenstates of the full Hamiltonian (1) are tensor products of the vibrational wavefunctions and the 
symmetrized qubit states. For bosonic atoms, permissible eigenstates are tensor products of external 
(i.e., vibrational) and internal (i.e., qubit) states of the same symmetry. Thus, anti symmetrized spatial 
wavefunctions are permitted for a pair of composite bosons, so long as their internal structure is also 
antisymmetric. 

We solve the Hamitonian (1) numerically for individual internal states. Examples of two-atom en- 
ergy spectra as a function of well separation are plotted in Fig. 2 for positive and negative interaction 
strengths a^j. We define E^^^ {d) to be the energy of the two-atom vibrational eigenstate with two- 

qubit internal state \(p) at well separation d. For example, the energy of |?/^^?/^^) |11) is El^-^^{d). 
As we have already discussed, not all combinations of n and \(p) are possible. For example, (d) 
is forbidden for identical bosons. This reduction of the size of the Hilbert space makes it possible to 
perform quantum gates adiabatically without losing coherence due to energetic degeneracies. Further- 
more, the Hamiltonian 's inherent particle and parity symmetries lead to selection rules which further 
enhance the fidelity of two-qubit operations. 

3. Universal entangling gate 

As shown by Hayes et al [11], it is possible to exploit these symmetries in order to produce a 
two-qubit entangling operation. We begin with a pair of identical atoms localized to opposite wells of 
the double- well Hamiltonian (1) in the far separated case (d a). We prepare the qubit on the left in 
the state |(/?") = (a|0) + /3|1)) and the qubit on the right in the state \(p^) = {fi\0) ^ . The initial 
wavefunction lipi), written as a tensor product of external and internal states, is then, 

1^,) = + 1^^) 1(^^(^")). (5) 
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Using 1^^) = ;^(|01) =t we rewrite equation (5) to make the symmetrization explicit: 

+ (|V^)+|V^))^(^|oo) + ^^|^+) + 5||ii)). 



(6) 



V2^ 



As the wells are brought together and separated adiabatically, the external states evolve according to 
Fig. 2. That is, each vibrational eigenstate at ^ a evolves continuously into its respective eigenstate 
aid = 0. As d decreases, the degeneracies between symmetric and antisymmetric eigenstates are lifted, 

resulting in a dynamic phase difference between + lip^ip^)) and -^(\^^^^) — \ip^ip^)), 

corresponding to the difference in respective energy curves (see energy curves 3e and 3f in Fig. 2a, or 
4c and 4d in Fig. 2b). Furthermore, degeneracies between the even two-qubit states {|00), |11), |^^)} 
are removed if the interaction strengths aij differ, which is usually the case. This state-dependent 
interaction results in additional phase differences between qubit states of the same symmetry [22]. 
Thus, each joint internal and external state acquires a unique phase, and the final state \ipf) upon re- 
separating the wells is. 



(7) 



'V2 

For positive scattering lengths, the phases are given by. 



1 r^' 



E^^..^ {d(t)) dt and 



1 f^' 



^,^4 m) dt. 



(8) 



Equation (7) is also valid for negative scattering lengths, although different phases will be acquired 
since ^(|^^^^) + |?/^^?/^^))is not the vibrational ground state when < 0. 

Clearly, this evolution can be thought of as the identity acting on the vibrational subsystem tensored 
with a unitary U acting on the qubit subsystem. Thus, we can discard the vibrational terms and examine 
the unitary evolution of the qubit subsystem by itself. Using matrix notation. 



|0) 



and 



(9) 



we can write U as, 



/2e~^^°° 
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(10) 
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Vo 

This entangling operation is diagonal in the partial Bell basis {|00), |^^), As noted in [11], 
even if the interaction strengths are state-independent, the singlet state |^~) acquires a phase differ- 
ent from the triplet states (except in the limit as aij ±oo, when the gate is no longer feasible). 
State-dependent traps or atomic interactions generally result in state-dependent interactions with the 
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environment as well, making the qubits more sensitive to noise. Avoiding this state-dependence leads 
to the inherent robustness observed in initial experiments [14], as compared to earlier experiments 
wherein gate fidelities were severely limited because of dephasing due to state-dependent traps [5]. 
Furthermore, this gate works for a wide range of positive scattering lengths, as we show elsewhere [22]. 
This is especially important for experiments employing atomic species with unknown or approximately 
known scattering lengths (e.g., Yb or Sr). 

A "controlled phase" gate (i.e., e~*^l^^^^^^l) in the computational basis can be obtained by com- 
bining single qubit phase gates S{0) = e~*^'^^^^' with a pair of U gates: 



G = U{S{7r)^S{0))U ■ 



\ 



-2i(f)Q, 









g-i(<t>++<t>-) 









-g-i(<t>++<t>-) 





(11) 
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G is locally equivalent to the "tunable controUed-phase" gate e , subject to the constraint: 

000 + 011 - 0+ - 0- = (2n± I) 7, V n G Z. (12) 

While a simple controlled-phase gate is itself a universal entangling gate, many quantum algorithms 
(e.g., the quantum Fourier transform) can be performed more efficiently when tunable controlled-phase 
gates are available. Here, the value of 7 can be easily tuned, simply by adjusting the speed at which the 
optical tweezers are combined and separated. The inherent robustness and easy tunability of this gate 
make it a highly desirable one for quantum information processing. 



4. Alternative entanglement preparation 

For some quantum information processing applications, universal entangling gates are not neces- 
sary and an ability to prepare entangled pairs will suffice. Atomic quantum repeaters based on entan- 
glement swapping [23] provide an example of one such application. We next examine a scheme which 
uses symmetrization requirements and a selective excitation to produce Bell pairs. 

We begin with two bosonic atoms in the ground state of a single well, both with internal state |0). 
The energy of this state is £'|^qq^(0). It is then possible to perform a coherent transition, selectively 
exciting to the eigenstate with energy ^|^^+^ (0). If the interaction for atomic qubits in the state |^+) is 
significantly different than that of |11), one can deterministically excite only to the |^+) qubit eigen- 
state, since the state 1 11) is off resonant and the overall state must remain symmetric. An excitation to 
the antisymmetric state |^~) is not possible as long as the two atoms are in the symmetric vibrational 
ground state. Thus, any initial population in antisynmietric vibrational states (e.g., due to heating) must 
be avoided to keep the fidelity of entanglement generation high. Furthermore, the vibrational spacing 
and sidebands due to the interaction energy must be spectroscopically resolvable. With typical vibra- 
tional energies on the order of kHz and on-site interaction shifts of close to lOOHz (dependent on 
atomic species but tunable by the tightness of traps), this selective excitation process is generally slow, 
but nevertheless viable. The final state after separating the atoms adiabatically is, 

l^final) = |(|^^) + |^^)) (|01) + |10)) . (13) 

This operation provides a novel way of creating Bell states deterministically, but does not constitute 
a universal two-qubit entangling gate. It does however allow for fundamental tests of quantum mechan- 
ics and Bell inequality violations, as well as basic quantum information processing and communication 
tasks. 
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Fig. 3. Snapshots of the magnitude of the two-atom vibrational wavefunction (||V^)| = | |V^(xa, Xb, t)) |) as a 

function of the position of each atom. Plot (a) shows the initial wavefunction, iV^init) = + \ip^i/j^))/ 

Plots (b-d) show the wavefunction after the wells have been brought together and separated. Initial conditions are 
the same for all figures, and only the well speed (in units of Vo — hcnJ^i^/Vo) is varied. The resulting vibrational 
state fidelities / = | (V^init|^) |^ are as follows: (b) v ^ O.Olvo, f = 0.9997. (c) v ^ O.lvo, f = 0.491. (d) 
v;^Vo,f = 0.002. 



5. Speed constraints for adiabaticity 

An approximation to the general adiabaticity criterion is given in [24] : 

« hLol V \a) + \h) , (14) 

where ujab = inin(|£^5(t) — Ea{t) \ /K) and where \a) and \h) are time-dependent eigenstates of an 
arbitrary Hamiltonian. Since our specific Hamiltonian is invariant under exchanges of both symmetry 
and parity, transitions between vibrational states of different symmetry or parity are suppressed. Thus, 
in our case ujah is determined by the energy gap of the two closest states having both equal symmetry 
and parity. This restriction contributes significantly to the robustness of this gate. Since only the double- 
well potential is time-dependent, the left side of the equation reduces to | ^^(^'^^(^)) | ^ Assuming constant 

V and maximizing | | ^ith respect to x, we obtain the adiabaticity criterion: 

V < ficjuolJVo. (15) 
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Time-dependent numerical simulations confirm the validity of this simple criterion over a wide 
range of values of Vq and aij (including a^j < 0). Plots of the two-particle wavefunction comparing 
adiabatic and non-adiabatic evolutions are shown in Fig. 3. Under adiabatic conditions, we recover both 
atoms in separate wells. Under non-adiabatic conditions both atoms may end up in the same well with 
non-negligible probability, resulting in an erroneous gate. However, time-dependent simulations have 
also shown significant revivals, with both atoms ending up in opposite wells with large probability even 
under non-adiabatic conditions. This suggests the very real possibility of producing a fast, coherent, 
non-adiabatic gate via optimal control. 

6. Conclusion 

In summary, we have proposed two schemes for preparing pairs of entangled atoms. We have shown 
it possible to construct a tunable universal entangling gate via the exchange interaction between iden- 
tical bosons, promising high fidelity operation for positive (repulsive) and even negative (attractive) 
interaction strengths. This is of particular importance for quantum information processing applications 
that use novel species of atoms, such at Group Il-like atoms (e.g., Yb and Sr) [22], for which the 
coUisional interaction parameters are partially unknown. In addition, we have introduced a novel en- 
tanglement scheme allowing the creation of Bell pairs. This scheme could prove useful for quantum 
communication schemes and fundamental tests of quantum mechanics. The use of this entanglement 
operation for Group Il-like atoms and its application to fundamental tests of quantum mechanics are 
studied in detail in other work [22]. 

7. Acknowledgments 

We specially thank D. Hayes and I. Deutsch for insightful discussions on entangling atoms via 
exchange interactions. We further thank Michael Skotiniotis and Nathan Wiebe for helpful comments 
and discussion. This work was supported by NSERC, AIF, CIFAR, iCORE, MITACS, NSF, and The 
Welch Foundation. 

References 

1. H.-J. Briegel, W. Diir, J. 1. Cirac, and P. ZoUer, Phys. Rev. Lett. 81, 5932 (1998). 

2. D. Gottesman and I. L. Chuang, Nature 402, 390 (1999). 

3. E. Knill, R. Laflamme, and G. J. Milburn, Nature 409, 46 (2001). 

4. R. Raussendorf and H.-J. Briegel, Phys. Rev. Lett. 86, 5188 (2001). 

5. O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. Hansch, and I. Bloch, Nature (London) 425, 937 
(2003). 

6. D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and R ZoUer, Phys. Rev. Lett. 82, 1975 (1999). 

7. E. Charron, E. Tiesinga, F. Mies, and C. Williams, Phys. Rev. Lett. 88, 077901 (2002). 

8. K. Eckert, J. Mompart, X. X. Yi, J. SchUemann, D. BruB, G. Birkl, and M. Lewenstein, Phys. Rev. A 
66, 042317 (2002). 

9. G. K. Brennen, I. H. Deutsch, and C. J. Williams, Phys. Rev. A 65, 022313 (2002). 

10. D. Jaksch, J. I. Cirac, R ZoUer, S. L. Rolston, R. Cote, and M. D. Lukin, Phys. Rev. Lett. 85, 2208 
(2000). 

11. D. Hayes, R S. Julienne, and I. Deutsch, Phys. Rev. Lett. 98, 070501 (2007). 

12. M. Kitagawa, K. Enomoto, KentaroKasa, Y. Takahashi, R. Ciurylo, P. Naidon, and P. S. Julienne, 
arXiv:0708.0752 (2007). 

13. M. Yasuda, T. Kishimoto, M. Takamoto, and H. Katori, Phys. Rev. A 73, 011403(R) (2006). 

14. M. Anderlini, P. J. Lee, B. L. Brown, J. Sebby-Strabley, W. D. Phillips, and J. V. Porto, Nature (London) 
448, 452 (2007). 

15. A. M. Dudarev, M. G. Raizen, and Q. Niu, Phys. Rev. Lett. 98, 063001 (2007). 

16. C.-S. Chuu, R Shreck, T. R Meyrath, J. L. Hanssen, G. N. Price, and M. G. Raizen, Phys. Rev. Lett. 95, 
260403 (2005). 



©2008 NRG Canada 



Babcock, Stock, Raizen, and Sanders 



9 



17. R. B. Diener, B. Wu, M. G. Raizen, and Q. Niu, Phys. Rev. Lett. 89, 070401 (2002). 

18. T. L. Gustavson, A. P. Chikkatur, A. Leanhardt, A. Grlitz, S. Gupta, D. E. Pritchard, and W. Ketterle, 
Phys. Rev. Lett. 88, 020401 (2007). 

19. J. Beugnon, C. Tuchendler, H. Marion, A. Gaetan, Y. Miroshnychenko, Y. R. P. Sortais, A. M. Lance, 
M. P. A. Jones, G. Messin, A. Browaeys, and P. Grangier, Nature Physics 3, 696 (2007). 

20. R. Stock, I. H. Deutsch, and E. L. Bolda, Phys. Rev. Lett. 91, 183201 (2003). 

21. T. Calarco, E. A. Hinds, D. Jaksch, J. Schmiedmayer, J. L Cirac, and R ZoUer, Phys. Rev. A 61, 022304 
(2000). 

22. R. Stock, N. S. Babcock, M. G. Raizen, and B. C. Sanders, arxiv: 07 11.0044 (2007). 

23. L.-M. Duan, M. D. Lukin, J. L Cirac, and R Zoller, Nature 414, 413 (2001). 

24. B. H. Bransden and C. J. Jochain, Quantum Mechanics, 2nd ed. (Pearson Education Ltd., Dorchester, 
Dorset, 2000). 



(c)2008 NRC Canada 



